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1.  INTRODUCTION 


1. 1  AIR  FORCE  GOALS  AND  NEEDS 

The  U.S.  Air  Force  needs  Improved  communications  methods  not  only  for  Increased 
bandwidth  but  also  for  Increased  security  and  margin  against  jamming,  fading, 
multipath,  and  Doppler  degradation.  The  purpose  of  this  study  was  to  determine 
which  tasks  In  the  receiver  operations  could  best  be  done  optically  and  to  see  how  to 
do  them.  The  particular  communication  method  we  studied  In  detail  Is  j»pread  spec¬ 
trum  communication  (hereafter  SSC)  and.  In  particular,  the  most  popular  SSC 
method  Involving  gseudo-random  codes  (PRC's)was  emphasized. 

There  are  many  possible  reasons  for  preferring  optical  signal  processing  to  digital 
signal  processing.  They  are  seldom  all  simultaneously  applicable.  Such  reasons 
Include 


•  parallel  processing, 

•  ability  to  handle  very  long  codes, 

•  high  bandwidth,  and 

•  compatibility  with  optical  storage  methods. 

We  believed  that  optical  processing  would  be  most  applicable  In  situations  where 
the  received  signal  Is  Itself  optical  (in  a  broad  sense)  and  noisy  conversions  of  sig¬ 
nal  (and  noise)  from  the  optical  domain  to  another  (electrical  or  different-optical) 
domain  could  be  avoided  during  the  first  stages  of  the  decoding  process .  We 
called  this  "direct  optical  decoding". 

1.2  ANALYSIS  OF  SPREAD  SPECTRUM  COMMUNICATION 


The  SSC  technique  utilizes  a  bandwidth  B.  to  transmit  a  message  of  bandwidth  B 
and  requires  m 


B.  >  >  B 


m 


(1-1) 


This  spread  In  bandwidth  (or  "spectrum")  offers  error  correcting  possibilities. 
Let  us  call 

R  «  Bt/Bm  (1.2) 

the  "spreading  ratio".  The  PRC  method  starts  with  a  string  of  R-blt*  code  words 
which  we  can  call  W-,  W„,  ....  W  .  With  this  string  of  code  words  we  can  send 
an  N-blt  message,  we  send  1,1..,  1  by  sending  W. ,  W„,  ....  W  .  We  send 
0,0 . 0  by  sending  Wj,  Wj . WN°,  where  1  N 

Wk  +  Wk°  =  0  (1-3) 

*These  bits  are  sometimes  called  "chips,  "  so  "bit"  can  be  used  for  our  "word". 
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for  k  =  1,  2,  ....  N.  Here  the  superscript  c  indicates  the  complement 

<1  *  -1,  -l6  =  1). 

There  are  several  possible  decoding  problems  with  the  PRC  version  of  SSC.  First, 
synchronization  to  better  than  Bf  *  is  required  between  the  received  signal  and  the 
code  string  being  used  for  decoding.  Even  when  synchronization  is  established,  it 
Is  not  automatically  maintained  because  transmitter- receiver  range,  atmospheric 
effects,  Doppler  effects,  etc.  affect  the  arrival  rate  of  signals.  Some  people  have 
suggested  that  extremely  long  codes  offer  an  extra  layer  of  security  in  SSC .  The 
time-bandwidth  product  for  the  total  code  word  string  is  RN.  Synchronization  with 
a  very  large  RN  product  might  discourage  enemy  attempts  to  decode  intercepted 
SSC  signals,  hi  addition  even  fairly  short  codes  can  offer  considerable  processing 
gain. 

Second,  even  with  synchronization,  decoding  may  be  difficult  because  of  atmos¬ 
pheric  effects  (multipath,  fading,  etc.)  and  operating  conditions  (Doppler,  poor 
tracking,  etc.).  These  problems  can  be  devastating  in  a  system  designed  to  be 
just  sufficient  in  their  absence.  Of  course,  sufficient  encoding  (including  burst 
error  correction)  can  work  in  severe  cases  but  the  bandwidth  suffers  even  more 
this  way.  Good  system  design  is  always  a  tradeoff. 

1.3  REPORT  OUTLINE 

We  deal  primarily  with  decoding  methods  (Section  2)  and  secondarily  with  a  new  SSC 
method  (Section  3).  Thereafter,  we  seek  to  assess  the  Impact  of  this  program 
(Section  4). 


2.  OPTICAL  DECODING  METHODS 

2 . 1  OPERATIONS  OPTICS  CAN  PERFORM 

We  distinguish  three  possible  tasks  for  optical  decoding: 

1  synchronization  with  long  codes, 

2  direct  optical  decoding  of  PRC  signals,  and 

3  error  correcting  decoding. 

2.2  SYNCHRONIZATION  WITH  LONG  CODES 

The  basic  concept,  explored  by  other  workers  earlier  (1),  is  to  use  the  two-dimen¬ 
sional  nature  of  the  optical  spatial  filter  correlator  to  correlate  the  received 
sequence  with  long  codes.  Two  problems  seem  to  appear.  First,  for  codes  of  the 
length  contemplated  just  storing  the  code  on  a  single  optical  filter  was  beyond  the 
state  of  the  art.  Second,  there  appears  to  be  a  bigger  problem  than  just  storing  a 
record  of  the  RN  code  bits  properly.  The  problems  which  affect  PRC  decoding 
over  the  time  of  one  code  word  affect  it  even  more  over  N  code  words.  Consider 
Doppler  effects,  for  Instance.  These  Introduce  subtle  temporal  scale  changes  in 
the  received  signal.  Furthermore  the  Doppler  effect  Is  neither  predictable  nor 
stable.  For  which  Doppler  shift  should  the  recognition  search  occur?  It  is  clear 
that  we  must  search  both  Doppler  shifts  and  codes  '  is  the  storage  problem 


b 


becomes  enormous. 


Let  us  accept  the  hypothesis  that  we  need  to  search  only  one  (probably  zero)  Doppler 
shift  and  draw  some  conclusions  and  determine  the  feasibility  from  the  storage- 
density  point  of  view. 

Horrlgan  and  Stoner  (1)  have  suggested  optical  correlation  for  code  synchroniza¬ 
tion  based  on  methods  derived  from  Cohen  and  Moses  (2).  This  has  the  advantage 
of  using  the  two-dimensions  available  to  optics  In  an  optimum  way.  Even  with  two 
dimensions  to  work  with,  the  total  code  length  Is  so  long  as  to  be  well  beyond 
current  technology. 

Let  us  imagine  a  code  of  several  days  length.  Periodically  the  sender  transmits 
the  message  l»l*U«*«l  (a  synchronizing  string).  Our  job  Is  to  locate  that  string  of 
M  bits  In  our  stored  record  of  P  bits,  where 

P  >>  M  >>  1.  (2.1) 

We  define 

Bt  =  base  code  transmission  bandwidth,  blts/sec.; 

T  =  M/Bt  =  duration  of  the  synchronizing  string,  sec.; 
t  =  l/Bj  *  duration  of  one  bit  or  chip,  sec.; 

T  =  P  /  B  t  =  duration  of  total  code ,  sec . ; 

-5 

D  =  t  x  1.157  x  10  =  duration  of  total  code,  days. 

We  assume  we  know  what  day  of  the  code  we  are  In.  We  divide  the  total  code  as 
shown  In  Figure  2.1.  We  should  have  no  doubt  as  to  which  one-day-long  segment 
It  belongs  to.  In  a  one-day  record  there  are 

Pt  =  8.64  x  104  Bt 


bits. 

Our  task  Is  to  utilize  the  methods  already  developed  for  optical  processing  and 
going  beyond  them  to  extend  their  applicability.  Thus  Horrlgan  and  Stoner  show 
that  handling  the  very  large  P  ^  values  required  for  one-day  codes  by  spatial  exten¬ 
sion  alone  Is  impractical  (see  their  Table  1). 
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1  A  four-day  code  (a)  becomes  a  cyclic  code  (b)  which  we 
represent  in  8  overlapping  one-day  sections  (c) 
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If  we  assume  that  21.5  cm  (4.6  cm  x  4.6  cm)  is  a  reasonable  mask  size,  we  can 
handle  0.36  minutes  at  B  =  10®.  This  is  a  factor  of  6000  off  from  the  D  =  lwe 
want.  The  question  then  Decomes:  "How  can  we  gain  three  orders  of  magnitude 
In  storage  without  increasing  the  device  size  and  complexity?"  Were  we  to  solve 
this  problem,  optics  would  offer  an  attractively  compact  and  efficient  vehicle  for 
solving  the  synchronization  problem. 

There  appears  to  be  at  least  one  way  to  accomplish  this:  use  wavelength  multiplex¬ 
ing.  Castro  et  al  (3)  have  developed  a  recording  material  which  allows  us  to  write 
over  1000  Independent  (orthogonal)  spectral  patterns  into  the  same  material.  A 
reprint  describing  their  method  is  attached.  Using  such  a  material  we  can 

(1)  process  the  Input  against  1000  masks  in  parallel  using  Incoherent  optical 
methods  and  read  out  the  detector  plane  in  sequence  using  heterodyne 
detection  or 

(2)  process  and  read  out  one  wavelength  at  a  time  using  either  coherent  or 
Incoherent  methods. 

As  both  processes  require  sequential  readout  and  thus  require  the  same  time,  we 
see  no  reason  to  do  parallel  processing  of  all  spectral  channels. 

The  question  now  is:  "Do  existing  lasers  give  sufficient  power  to  do  this?" 

During  T  (assumed  for  the  moment  to  be  0.36  minutes),  we  may  want  to  process 
1000  "frames"  In  separate  spectral  channels.  We  then  have 

T  =  T/1000  (2.2) 

r 

as  the  time  to  Interrogate  each  frame.  For  our  example,  Tp  =  22  msec.  Follow¬ 
ing  Horrigan  and  Stoner's  example,  we  have 

M  =  2.147  x  109. 

Even  if  we  use  a  CCD  type  detector  array  (none  this  size  exist),  we  cannot  Interro¬ 
gate  at  the  rate  B  much  less  at  MB.  What  we  need  Is  a  way  to  reduce  the  electronic 
processing  rate  substantially.  The  only  way  we  see  to  do  this  is  to  use  an  optical 
bistable  device  (4)  to  threshold  out  all  but  the  correlation  due  to  perfect  synchroni¬ 
zation.  We  would  have 

tp  =  0.36  x  10  3  min.  =  19.6  msec 

to  determine  the  x-y  location  of  that  peak  when  It  occurs.  Centroid  detectors 
with  greater  than  Vm  =  resolvable  units  In  each  direction  are  available  commer¬ 
cially  (5).  These  are  silicon  detectors,  so  we  can  assume  a  good  quantum  effici¬ 
ency,  t) .  We  will  assume,  pessimistically,  ^  =  0.1.  To  locate  to  one  part  in  Vm\ 
we  want  at  least 


N  =  100  M 
P 

photons.  Assume  a  wavelength  \ ,  and  a  laser  power  PT  then  we  want 


>  100  M. 

where  h  is  Planck's  constant  and  c  Is  the  speed  of  light.  For  \  =  0.5/im  (green 
light) ,  we  need  (for  this  example) 

PL>  100  x  2.147  x  109  x  6.626  x  10~34  x  3  x  IQ8 
0.1  x  0.5  x  10-6  x  1.96  x  10-2 

=  43.5  |jw. 


This  level  is  easily  obtained  with  every  commercially  available  tunable  dye  laser. 
The  frequency  tuning  rate  must  be  faster  than  20  msec  per  wavelength.  One  milli¬ 
second  times  are  not  unreasonable  for  a  stepping  motor  to  rotate  a  prism  to  change 
the  wavelength. 

It  thus  appears  that  the  combination  of  processing  techniques  already  devised  with 
some  of  the  newest  optical  devices  (wavelength  hole  burning,  bistable  devices,  and 
centroid  locators)  can  handle  the  full  synchronization  problem.  This  system  is  very 
compact  and  can  perform  correlations  over  the  whole  one-day  code  during  a  time 
equal  to  that  required  to  accumulate  M  bits,  so  all  bits  can  be  tested  on  line. 

A  far  more  modest  task  can  be  valuable  as  well.  Most  digital  electronic  correla¬ 
tors  can  handle  no  more  than  about  1000  bits,  while  optics  can  handle  10,  000  or 
more  rather  easily.  Decoding  "short"  codes  may  represent  a  useful  application 
of  optics. 

2.3  DIRECT  SPECTRAL  DECODING 
2.3.1  Overview 

We  have  reached  here  the  key  part  of  our  effort.  The  idea  is  to  postpone  optical  - 
to-electrlcal  conversion  until  after  PRC  decoding  (with  various  sources  of  noise 
and  jamming  filtered  out  optically  and  Doppler,  synchronization,  fading,  and  point¬ 
ing  errors  continually  counteracted  optically).  Thus  the  detector  sees  only  the 
noise-filtered  message  at  bandwidth  Bm-  Large  gains  in  jam  immunity  and  in 
signal -to-noise  ratio  may  be  expected  If  such  a  task  can  be  accomplished. 

A  question  naturally  arises  as  to  what  wavelength  limitations  exist  in  this  method. 
We  shall  see  that  the  limitations  are  largely  on  materials  in  the  components  and, 
in  particular,  in  acousto-optic  materials  which  we  know  to  work  from  0.3  pm  to 
11  pm.  We  have  no  data  on  longer  wavelengths,  so  the  ultimate  long  wavelength 
limits  are  unknown. 

Figure  2.2  shows  the  major  system  components.  A  standard  tracking  telescope 
keeps  the  centroid  of  the  signal  through  the  narrow  band*  filter.  This  beam  is 
then  spectrally  filtered  by  a  new  acousto  optic  tunable  niter  (AOTF's)  Invented 


*  Wide  enough  to  encompass  all  reasonable  Doppler  shifts. 
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under  this  contract  to  achieve  a  narrower  spectral  band  pass  than  prior  AOTF's 
and  to  track  Doppler  shifts.  The  beam  so  tracked  and  filtered  passes  through  a 
t_lme  Integrating  correlator  (TIC)  to  a  detector  array  which  allows  synchronization 
tracking.  Fading  (measured  after  spectral  filterlng)and  Doppler  shifting  (measured 
by  the  correction  needed  to  maintain  the  proper  location  of  the  spectrally  filtered 
signal)  are  taken  into  account  in  the  amplitude  and  time  scale  of  the  signal  fed  to 
the  TIC. 

Finally,  the  signal  correlated  against  is  not  a  section  of  the  code  word  string 
Wi>  W2,  ....  but  a  modified  version  far  more  tolerant  of  Doppler,  etc. 

This  modified  correlation,  which  we  call  "generalized  matched  Hlterlng"  (GMF) 
is  certainly  the  most  wldely-appllcable  invention  of  this  work. 

2.3.2  The  Super  High  Resolution  AOTF 

Here  we  discuss  the  AOTF  system  Invented  under  this  contract.  We  believe  it  is 
widely  applicable.  AOTF's  have  become  valuable  tools  in  modern  electrooptic 
systems  (6).  This  invention  is  a  new  way  to  use  them. 

In  their  simplest  form,  AOTF's  are  simply  variable-spacing  gratings.  In  fact  they 
are  thick  holographic  gratings,  so  their  Bragg  diffraction  efficiency  can  be  quite 
high  (7).  One  difficulty  with  AOTF's  has  been  that  their  resolving  power  is  still 
rather  modest.  We  show  here  a  rather  simple  way  to  Increase  their  resolving 
power  significantly  at  a  concomitant  price  of  reducing  the  tunabllity  range.  Indeed 
the  number  of  spectral  channels  is  conserved  but  their  width  Is  now  controllable 

The  basic  idea  is  to  achieve  the  high  dispersion  with  a  nontunable  element,  e.g. 
a  grating.  The  tuning  is  done  by  changing  the  angle  of  incidence  of  light  on  the 
fixed  disperser  with  an  AOTF.  Together,  the  AOTF  and  the  fixed  disperser  com¬ 
prise  a  system  equivalent  to  a  single  highiy-dispersive  grating  with  limited  tunability. 

In  the  analysis  which  follows  we  use  simple,  thin-grating  equations  for  simplicity. 
Figure  2.3  shows  a  simple  grating  of  length  L  and  spacing  d  diffracting  light  of 
wavelength  X  from  an  incident  angle  8  to  an  exit  angle  <p  .  For  constructive  inter¬ 
ference  we  require 

d  sln8  +  d  sln0=  nA,  (2.4) 

where  n  is  an  integer  (the  diffraction  order) .  In  other  words 

sin0  =  _  sln0  (2.5) 

Diffraction  limits  the  resolution  to 

Asln0  ~  .  (2.6) 

The  effect  of  dispersion  is 

Asin0  =  .  (2.7) 
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Thus  the  diffraction-limited  resolving  power  is 

«  il)  *  "  V 

where 

NQ  =  L/d 


(2.8) 


(2.9) 


is  the  number  of  lines  in  the  grating.  Figure  2.4  shows  two  gratings  of  length 
L  cascaded.  If  the  grating  spacing  of  the  first  grating  is  d  and  the  grating  spacing 
of  the  second  grating  is  S  ,  it  is  easy  to  see  that  the  output  angle  \f/  is  given  by 

sin  tjf  =  ^7-  +  sin  0  .  (2.10) 

Here 

1_  1  1  . 

S'  “  S  ~  d  (2.11) 


Calling  «M\,d')  the  angle  we  would  observe  for  wavelength  a  if  8  =  S',  we 
have 


sin  U,  S')  -  slni|»(x,S)  -  -j-  . 

We  can  tune  d  over 

L/N  <:  d  c  L. 

o 

This  swings  sin  if/  (A .  S ')  over 

YT  *  sin  ^  (X. 8)  -  sin  0  (A, 6')^ 

Dispersion  gives 

A  sin  (X,  S')  *  A  sin  (a, 8)  = 
so  the  wavelength  tuning  range  is 


(2.12) 


(2.13) 


(2.14) 

(2.15) 


For  practical  cases  we  might  have  d/g  ~  100,  so  we  can  get  a  factor  of  100 
increase  in  resolution  at  the  price  of  a  factor  of  100  decrease  in  the  tuning  range. 

It  Is  obvious  that  such  an  AOTF  system  could  be  very  useful  in  decoding  a  frequency 
hopping  system  with  only  one  channel  rather  than  one  channel  per  frequency. 
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Using  this  new  AOTF  system  we  seek  to  keep  the  maximum  signal  on  axis.  From 
the  tuning  required  to  do  so  we  track  Doppler  shifts .  We  Integrate  over  a  time  of 
about  N/Bt  to  assure  a  sample  Independent  of  the  exact  form  of  W^. 

2.3.3  Time  Integrating  Correlator 

The  "classical"  TIC  (8)  (9)contemplates  modulating  a  laser  beam  to  create  a  time 
signal  Si  (t),  superimposing  Si  (t)  at  all  points  along  an  acousto-optic  modulator 
(x  direction  along  Its  axis),  sending  the  signal  S2  (t  -  *)  across  the  modulator  at 
speed  V,  Imaging  the  modulator  onto  a  detector  array,  and  time  Integrating.  The 
output  signal  at  x  Is 


0  (x. 


ll’ 


S1(t)S2  (t  -  5)  dt. 


(2.17) 


which  we  recognize  as  a  cross  correlation  function  If  Si  and  S2  are  real. 

In  our  case  we  already  have  a  modulated,  monochromatic  light  beam,  so  all  we 
need  to  do  Is  Impact  It  on  an  acousto-optic  modulator  driven  by  S2  (t)  and  Image 
onto  the  detector  array.  We  should  obtain  the  peak  signal  (Integrated  over  a  time 
of  N/Bt)  at  the  center  of  the  array  for  perfect  synchronization.  Now  we  can  track 
synchronization  continuously  with  a  time  constant  of  N/Bt- 

2.3.4  Generalized  Matched  Filters 

Here  we  reach  what  we  believe  to  be  the  most  Important  invention  resulting  from 
this  work.  In  the  following  description  we  cast  the  GMF  In  terms  of  Its  two-dimen¬ 
sional  aspect  as  well  as  Its  one-dlmenslonal  aspect.  We  conclude  that  wherever 
matched  filters  are  useful,  generalized  matched  filters  are  probably  more  useful. 

Matched  filtering  has  become  the  primary  tool  of  optical  pattern  recognition 
despite  the  fact  that  matched  filters  fall  to  take  into  account  either  the  differences 
among  patterns  of  classes  we  wish  to  distinguish  or  differences  among  patterns 
belonging  to  the  same  class.  Thus  a  filter  matched  to  the  letter  "A"  Is  the  same 
whether  or  not  we  wish  to  tolerate  rotations,  scale  changes,  or  other  distortions 
and,  of  course.  Independently  of  the  magnitudes  or  probabilities  of  those  distortions. 
Likewise  the  filter  matched  to  the  letter  "A"  Is  the  same  regardless  of  what  we  are 
trying  to  distinguish  It  from.  e.g.  "B",  "C"  or  "B"  and  "C".  For  many  years 
we  have  sought  a  way  to  adapt  filters  to  the  specific  job.  Caulfield  and  Maloney  (10) 
showed  that  filters  matched  to  linear  combinations  of  the  undlstorted  "character 
alphabet"  could  be  made  to  give  maximum  dlstlnguishabillty  between  undlstorted 
Input  patterns.  They  also  showed  that  these  new  filters  were  somewhat  more 
tolerant  to  Input  distortions  than  were  matched  filters.  What  we  describe  here  Is 
a  fully  general  solution  to  the  problem  of  generating  the  optimum  filter  (which  of 
course,  may  not  be  the  matched  filter)  for  each  particular  problem.  We  call  these 
new  filters  "generalized  matched  Alters"  or  GMF's  because  they  reduce  to  the 
matched  filters  (MF's)  when  the  problem  Is  the  simple  one  for  which  MF's  are 
known  to  be  optimum.  For  all  more  complicated  cases  the  GMF's  will  be  superior 
to  MF's  In  both  Insensitivity  to  distortions  we  want  to  tolerate  and  sensitivity  to 
between-class  differences. 
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While  the  GMF  Is,  to  our  knowledge,  a  new  concept,  the  mathematical  tools  upon 
which  it  is  based  are  very  old.  The  new  contribution  described  here  is  that  of 
formulating  the  GMF  problem  in  such  a  way  that  it  is  clear  that  the  prior  mathema¬ 
tical  methods  describe  it. 

Let  us  describe  digital  matched  filtering.  A  discrete-point  Input  function  f(xm,yn) 
is  Fourier  transformed  to  produce  a  discrete-position  Fourier  transform  F(uk,  Vf  ). 
Then  F(uk,vf  )  is  multiplied  by  the  matched  filter  M(uk,  vf  ),  that  product  is 
retransformed,  and  the  resulting  pattern  g(5p,  rjq)  is  examined.  If  the  centroid 
of  f(xm,yn)  is  at  xm  =  0,  yn  =  0;  then  we  want  to  evaluate  f(0, 0)  is  simply 


g(0,0)  -  EE  F(uk,v.  )M(uk.v#  ).  (2.18) 

k  f 

Formulated  in  this  way,  the  problem  of  deriving  the  filter  M(uk,  v()  Can  be 
recognized  as  a  special  case  of  "linear  discriminant  analysis"  (H)whlch,  in  turn, 
is  a  special  case  of  "principal  component  analysis"  (12).  Rather  than  repeat  the 
known  mathematics  In  detail  here,  we  simply  describe  what  linear  discriminant 
functions  (LDF's)  are  and  in  what  sense  they  are  optimum.  We  suppose  that 
N  measurements  are  made  on  each  input  object  and  call  those  measurements  for 
the  unknown  object  x  the  measurement  vector  1c  =  (xj,  Xo,  . . . ,  Xvj).  Let  there  be 
K  classes  of  objects  (Cj,  ....  Cjr)  and  K  linear  discriminant  functions 
(LDFj,  . . . ,  LDFk)-  The  operation 

LDFj  (x)  =  S!(x)  (2.19 

gives  a  real  number  Sj  (x) .  Indeed 

LDFt(x)  =  Vj  *X  +  Soi,  (2.20) 

where 

Vi  =  (V* . VN)  (2.21) 

is  a  set  of  real  numbers  and  S0i  Is  a  real  number. 

The  LDF's  are  so  chosen  that 

[LDFj  (x  e  Cj)  -  LDFj  (x  rfCj)]  (2.22) 

is  maximum,  where  £  [•  ]  is  the  expected  value  operator.  LDFi  is  that  linear 
combination  of  x-'s  most  likely  to  distinguish  between  x  «  Ci  andx  4  C|. 

Usually  the  LDF's  are  normalized  so  that 

£  LDFj  (x  t  cj>  1=  8ij  (2.2S) 
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Thus,  for  an  unknown  x,  we  evaluate  all  of  the  S[  (x)'s.  If  the  maximum  one  is 
Sj[(x),  we  assign 

x  (  Ck.  (2.24) 

The  optimization  Is  designed  to  give  a  minimum  expected  error  using  such  a 
decision  algorithm. 

Our  approach  Is  to  make  the  feature  vector  for  f(xm,  yp)  the  set  of  {  FOty,  v()  } 
ordered  In  any  arbitrary  way.  The  LDFi's  give  a  "filter"  VO^,  Vj).  The 
S0  fa  give  thresholds. 

Let  us  now  look  at  some  GMF's.  If  we  expect  no  distortions  of  f(s,y)  and  no  other 
Images  In  the  field,  then  we  can  write 


Cx  =  |f(,y)  +  n(x,y) } 

(2.25) 

and 

C2  =  {n(x,y) } 

(2.26) 

where  {n(x,  yj)  represents  the  noise 

For  this  case  we  would  expect 

LDFj  =  MFt 

(2.27) 

and,  If  £[n(x,y)]  =  nQ, 

LDF2  =  -  MFj 

(2.28) 

That  Is,  to  recognize  noise  (C2)  we  say  "noise  is  what  Is  not  signal", 
simplicity  we  chose 

For 

C.  =  { rect(x)  +  n(x)} 

(2.29) 

and 


C2  =  { n(x)}  (2.30) 

where 


£[n(x)]  =0.1. 

Figure  2.5  shows  the  two  GMF's.  We  recognize 
GMF1  =  MFj  =  sine  (u) 


(2(31) 


(2.32) 
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and 


GMFg  =  MF2  *  -  sine  (u) 

Thus  the  MF  Is  a  special  case  of  the  GMF 
Now  If  we  have 
Cj  *  { rect(x)  +  n(x)} 

C2  *  { rect(2x)  +  n(x)}  ,  and 
Cg  =  {n(x,y)}  , 

we  have  no  expectation  other  than  that  the  GMF's  should  not  be  the  MF's. 
Figure  2.6  shows  the  resulting  GMF's. 

To  test  performance  of  the  GMF's  we  went  to  an  even  harder  problem: 


(2.33) 


(2.34) 

(2.35) 

(2.36) 


=  {  £  rect  (6x)  +  n(x)  }  , 

(2.37) 

C2  =  {  ~  rect  (7x)  +  n(x)  }  , 

(2.38) 

c3  =  {  |rect  (8x)  +  n(x)  }  . 

(2.39) 

C4  -  n(x). 

(2.40) 

We  placed  samples  of  these  noisy,  equal -energy  rectangle  functions  In  a  series 
as  shown  In  Figure  2.7.  Operation  on  this  series  with  the  first  three  MF's 

Figure  2.8  gave  poor  discrimination.  Operating  on  the  same  series  with  the 
first  three  GMF's  Figure  2.9  gave  excellent  discrimination.  Measuring  at  the 
centers  [  g(0,0)  of  our  previous  analysis]  we  obtained  the  following  response 
matrices  (normalized  to  unity  diagonal  elements): 

MF  RESPONSES 

1  2 

3 

N  1  /  1  096 

0.92  \ 

U  2  1  0.81  1 

J  3  \  0.64  0.82 

i  / 
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Figure  2.6  The  GMF’s  for  two  rectangle  functions  in  white  noise  and  white  noise. 
Note  that  slnc-Like  appearance  vanishes. 
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The  response  of  the  pulse  train  of  Figure  2,7 

operated  on  by  MF's  for  (a)  the  highest  (narrowest)  rectangle. 

H  (b)  the  middle  height  rectangle,  M,  and  (c)  the  lowest  (widest) 
rectangle  L.  Note  the  ambiguities,  especially  with  L. 


The  response  of  the  pulse  train  of  Figure  2.7  to  GMF's  for  (a)  H, 
(b)  M,  and  (c)  L.  Note  the  resolution  of  all  ambiguities  as  well  as 
the  "anticorrelations"  at  the  centers  of  the  "wrong  pulses". 
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GMF  RESPONSE 


The  average  value  of  the  difference  between  on-axls  and  off-axis  elements  is  only 
0.14  for  MF's  (poor  discrimination)  but  is  1.10  for  GMF's  (good  discrimination) 
when  the  on-axis  terms  are  equal.  Another  comparison  is  worst  cases.  The  worst 
case  differences  for  MF's  is  0.00.  The  worst  case  difference  for  GMF's  is  0.62. 

While  many  practical  problems  with  optical  implementation  of  GMF's  are  obvious 
(separating  real  and  imaginary  parts,  digital  filter  generation,  handling  large 
numbers  of  GMF  sample  points,  etc.),  the  basic  conclusions  are  clear.  Matched 
filtering  is  a  special  case  of  the  far  more  general  and  powerful  technique  which  we 
call  generalized  matched  filtering. 

2 . 4  ERROR  CORRECTING  DECODING 

While  the  method  we  describe  is  both  simple  and  promising, its  explanation  can  be 
confusing.  Careful  definition  of  the  problem  is  essential.  We  receive  MN  message 
bits,  where  M  and  N  are  preassigned  Integers.  The  bits  arrive  serially  and  can  be 
called  bi,  b2,  ....  bjyiN-  we  have  a  burst  error  correcting  code,  we  must 
rearrange  the  received  bits  before  decoding.  Typically,  we  might  arrange  them 
into  N  groups  of  M  bits  each  something  like: 


bl  bN+l 
b2  bN+2 


bNM+ 


1-N 


bN  b2N 


bNM 


Thus  a  burst  error  of  up  to  N  bits  can  affect  at  most  one  bit  in  each  of  the  N  words 
(rows).  Now  if  we  can  correct  for  a  few  (much  less  than  M)  bit  errors  in  a  word, 
we  have  achieved  a  large  degree  of  Immunity  to  burst  errors.  To  correct  for  up 
to  P  bit  errors,  we  restrict  the  transmitted  message  to  being  words  from  the  set 
of  codes  {  CW^  CW2,  . . .  CWq}  which  have  the  property  that  for  i  ft  ),  CWj  and 
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CWj  differ  In  at  least  2FM-1  places  (the  "Hamming  distance"  Is  at  least  2P+1).  Iden¬ 
tifying  each  row  of  the  matrix  just  formed  (e .  g .  row  1  comprising  bj,  .  .  • 
as  a  reoetved  word  RW^  (RW^  In  this  example),  we  want  to  find  the  Hamming 
distance  between  RWjj  and  each  of  the  Q  code  words .  The  closest  code  word  Is 
assumed  to  be  the  message  word.  Clearly  P  bit  errors  can  be  corrected  in  this  way. 

We  now  recognize  that  If  we  encode  l's  and  0rs  as  spatially-complementary  patterns, 
e.g. 

pfT  and  tM  , 

then  template  matching  (or,  equivalently,  matched  filtering)  of  the  received  word 
with  each  code  word  will  give  signals  whose  differences  are  proportional  to^the 
Hamming  distances  between  the  received  word  and  each  of  the  code  words.  Thus 
matched  filtering  can  be  used  for  error  correcting  decoding. 

The  required  process  is  shown  in  Figure  2.10,  while  Figure  2,11  shows  the  system 
layout.  In  Figure  2.10,  we  show  five  steps.  Steps  3  and  4  can  be  done  optically 
as  shown  in  Figure  2.11.  Each  received  word  is  arranged  on  one  line  in  the  x 
direction.  A  total  of  N  received  words  can  be  handled  in'parallel.  A  cylindrical 
lens  Fourier  transforms  X  into  U  and  allows  spreading  in  the  V  direction.  The 
holographic  mask  has  Q  multiplexed  words  encoded  on  it,  each  at  a  distinct  angle 
and  all  In  the  U  direction.  There  is  no  V  variation.  A  second  cylinderical-spherical 
lens  Fourier  transforms  the  U  direction  into  £  and  images  V  into  rj .  Along  the 
output  £  row  corresponding  to  the  input  x  rows  of  code  word  n ,  we  find  Q  correlation 
spots  (one  for  each  of  the  error  correcting  code  words.  The  £  value  corresponding 
to  the  strongest  correlation  gives  the  error  corrected  message  of  word  n. 

Optical  decoding  has  two  great  advantages  over  electronic  decoding  -  speed 
(from  parallel  operations)  and  versatility  (all  codes  are  equally  easy  to  handle. 

2 . 5  CONCLUSIONS  OF  THE  DECODING  STUDY 

It  13  now  clear  that  spread  spectrum  signal  decoding  by  optical  means  Is  practical 
and  useful  for 

.  synchronizing  with  longer-than-1000-blt  codes  (if  the  code  lengths  are  not 
so  long  as  to  preclude  synchronization  altogether), 

.  correlating  In  a  way  designed  to  take  statlstically-predtctable  distortions 
into  account,  and 

.  decoding  burst -error-correcting  codes . 

The  direct  optical  decoding  of  pseudo-randomly  coded  spread  spectrum  radiation 
(so  that  the  detectors  see  only  the  decoded  message  bits)  Is  both  feasible  and 
worthwhile . 

♦This  follows  directly  from  the  template-matching  Interpretation  of  matched  filtering. 
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STEP  I 


STEP  2 


STEP  3 


STEP  4 


Figure  2.10  Desired  Operations  For  Optical  Error 
Correcting  Decoding 
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2.11  Optical  Arrangement  For  Error  Correcting 


3. 


A  NEW  SPREAD  SPECTRUM  METHOD 


3.1  INTRODUCTION 

The  basic  Idea  here  is  to  use  polychromatic  radiation  which  is  spectrally  encoded. 
One  such  method  is  called  "frequency  hopping".  In  frequency  hopping,  successive 
message  bits  are  sent  over  discrete,  non-overlapping  spectral  channels.  Thus  the 
bandwidth  for  the  message  exceeded  that  for  any  bit.  Our  approach  was  to  use  the 
full  bandwidth  for  each  bit.  The  way  to  do  this  while  still  utilizing  spectral  encoding 
is  to  ullllze  different  spectral  weightings  for  different  "channels".  For  example, 
we  might  operate  on  the  Fourier  transform  of  the  source  spectrum  in  such  a  way 
that  the  various  "channels"  are  separable  only  in  Fourier  transformation. 

3.2  SPATIAL  FOURIER  TRANSFORMATION 

Here  we  want  to  show  how  a  spatial  display  of  the  cosine  Fourier  transform  of  the 
source  spectrum  can  be  obtained  by  interferometry.  The  basic  insight  is  due  to 
Michelson  and  Is  covered  well  by  Born  and  Wolf  (13) .  Modern  interferometers  to 
accomplish  this  are  reviewed  by  Caulfield  (14.  15,  16).  The  basic  idea  is  that  an 
interferometric  fringe  pattern  is  simply  the  spatial  display  of  modulo  X  (X  being 
the  wavelength)  equal -paths  from  the  source.  If  the  path  difference  is  a  linear 
function  of  position  in  the  fringe  field  (as  when  the  fringes  are  formed  by  two 
mutually-inclined  collimated  beams),  the  envelope  of  the  interference  pattern  is 
the  cosine  Fourier  transform  of  the  spectral  source  distribution.  Figure  3.1 
shows  an  example  from  Born  and  Wolf  (IS)  who  plot  the  envelope  or  fringe  visibility 
as  a  function  of  source  bandwidth  in  wavenumber  (which  they  call  Ak)  and  path 
difference  (which  they  call  A  V). 

Consider  an  Interferometer  which  splits  the  light  into  two  collimated  beams  and 
Interferes  them  at  an  angle.  Let  x  be  the  dimension  normal  to  the  bisector  of  the 
angle  between  the  two  beams  and  in  the  plane  of  that  angle.  Call  x  =  o  the  position 
of  perfect  path  matching.  At  all  other  points  (x  *  o)  the  wavefront  for  one  beam 
passes  x  before  the  wavefront  from  the  other  has  arrived.  The  time  difference 
is  proportional  to  x.  Thus  we  have  a  linear  spatial  display  of  relative  time 
differences  between  the  rays  travelling  the  two  paths.  The  power  between  wave- 
numbers  cr  and  a  +  do  is  S(o  )dcr  .  The  amplitude  of  that  light  varies  as  the  cosine 
of  the  frequency  v  -  ca  times  the  relative  time  delay  (proportional  to  x).  The 
sum  (Integral)  of  3(cr )  over  all  u 's  is  simply  the  cosine  Fourier  transform  of 
S(o). 

In  earlier  applications  (14)  this  spatial  pattern  was  used  as  a  means  for  deducing 
S(ct  )  rather  than  measuring  it  directly  with  a  dispersive  spectrometer.  Because 
no  scanning  occurs,  arbitrarily  short  pulses  can  be  handled  in  this  way.  Because 
all  wavenumbers  contribute  at  every  point,  some  dynamic  range  problems  are 
obviated.  The  application  we  will  describe  here  is  totally  new.  This  is  the  first 
study  devoted  to  spectral  modulation  in  such  an  Interferometer. 
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3.3  SPATIAL  FILTERING  OF  SPECTRAL  FOURIER  TRANSFORMS 


Let  us  consider  two  modified  Mach-Zehnder  Interferometers  In  series  which  have 
compensating  effects .  That  Is  all  paths  through  the  two  are  of  equal  length,  but 
there  Is  an  Intermediate  plane  where  the  cosine  Fourier  transform  appears. 
Suppose  that  the  source  spectrum  Is  S(cr),  Then  the  spatial  cosine  Fourier  trans¬ 
form  Is  „ 


s(x) 


=  J  S(<r) 


S(o  )  1  +  cos  (2mr  xcosd)  do  , 


(3.1) 


where  a  Is  the  wavenumber  and  S(o  )do  Is  the  power  between  o  and  o+do  . 
But  we  bound  the  x  plane  between  Xj  and  X2 ,  so  the  output  spectrum  is  (except 
for  Irrelevant  multiplicative  factors) 
x„ 


’•<'>  -/ 


/  8(X) 

1  +  COS  (2ir  or  xcos# ) 

dx. 

(3.2) 

X 

- 

If  we  Insert  a  mask  M(x)  In  the  Fourier  transform  plane  we  obtain  an  output 
spectrum 


Sm<ff> 
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[ 


s(x)  m(x)  1  +  cos(2iro  xcosfl) 


1  dx. 


J 


(3.3) 


The  cosine  Fourier  transform  of  that  spectrum  Is 


s  (x) 
m'  ' 


CD 

-/  V"  >  [ 


Sm^a  )  |  1  +  cos  (2ir  o  xcos#  ) 


do. 


(3.4) 


The  display  sm(x)  Is  the  result  of  Insertion  of  the  mask  m(x).  It  was  our  Initial 
hope  that  we  would  find 


sm(x>  *  B+m(x), 


(3.5) 


where  B  Is  a  background  term.  In  retrospect  that  hope  can  be  seen  to  be  naive. 
The  best  example  of  the  proper  approach  Is  the  Selective  Interferometric  Modu¬ 
lation  System  (SIMS)  which  Is  a  special  case  of  our  method  (17).  In  SIMS  we  set 


m(x)  *  1  +  cos  (2ir o  x  cob  6).  (3.6) 

Note  that  the  peaks  of  m(x)  hit  where  the  peaks  of  1  +  cos  (2ir  o0  x  cos  $ )  are. 


31 


The  effect  is  to  transmit  light  of  wavenumber  aQ  preferentially.  Then  they 
move  to 

m(x)  ■  •  1  +  sin  (2ir  cr  0  x  cos#)  (3 . 7) 

which  attenuates  selectively.  Thus  SIMS  can  modulate  one  wave  number 
component  of  the  spectrum.  We  must  modulate  in  a  well-defined,  distributed 
filter  function.  The  mask  which  leads  to  a  peak  at 

x  =  I/ct  cos  6  (3.8) 

o 


is,  of  course. 


m(x)  =  1  +  cos  (2-jt <to  x  sin  6  ).  (3.9) 

Thus  it  Is  clear  that  sm(x)  and  m(x)  are  not  similar  in  shape.  Furthermore  the 
example  we  chose  [  sm(x)  with  a  local  peak  at  x  =  l/aQ  cos  $  ]  is  easily  detectable 
with  a  spectrometer  tuned  to  cr  =  cr0.  A  communication  system  based  on  this  is 
no  more  secure  than  a  system  based  on  a  prism  or  grating  spectrometer.  What  we 
need  is  to  choose  a  more  complicated  x  pattern  we  want  to  detect.  Call  it  p(x). 
Then  we  set 


m(x)  = 


P(x) 


s(x)  +  C  ’ 

where  C  is  a  small  constant  used  to  avoid  singularities. 

Now  it  is  clear  that  we  can  choose  N  orthogonal  patterns  p1(x),  p  (x), 
That  is  1  1 

2 


1. 


Pj(x)  Pj(x)  dx  =  <5  ^ 


(3.10) 


...  PN(x). 


(3-11) 


The  corresponding  masks  are  mj(x),  m2(x),  ....  mn(x).  We  can  use  N  independent 
channels  separated  in  the  direction  normal  to  the  beam  splitting  plane.  By  modu¬ 
lating  the  k rty  channel  and  observing  through  the  corresponding  mask  PkWLi1}  the 
spatial  Fourier  transform  of  the  receiver,  we  can  communicate  over  the  ktft/ 
channel  independently  of  communication  over  any  other  channel.  Furthermore, 
without  knowing  P]c(x),  we  are  likely  to  detect  some  combination  of  all  N  messages. 
Thus  we  achieve 

»  multiplexing  of  N  signals, 

•  orthogonality  of  those  signals 

,  security  (the  would-be  code  breaker  must  decide  that  this  is  a  Fourier 
transform  code,  examine  the  proper  range  Xj  <:  x^  X2>  and  read 
through  the  proper  mask),  and 

v  jam  immunity  (jamming  with  either  narrow  band  or  broad  band  radiation 
adds  only  to  the  background  but  should  not  affect  the  signal  modulation. 
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The  main  drawbacks  are 


•  great  complication, 

•  need  for  bright,  broadband  radiation  sources,  and 

s  relatively  low  depth  of  modulation  and  hence  great  power  needs. 

3 . 4  EXPERIMENTAL  TESTS 

The  setup  we  wished  to  accomplish  Is  shown  In  Figure  3.2.  The  actual  setup  is 
shown  schematically  In  Figure  3.3.  The  upper  part  of  Figure  3 . 3  represents  the 
transmitter.  The  two  Interferometers  of  the  transmitter  are  combined  into  one 
Interferometer  traversed  twice  by  Insertion  of  the  reversing  mirror  Mx.  Clearly, 

If  Mx  Is  Inserted  properly,  all  rays  from  the  source  to  the  output  lens  Lz  have 
equal  path  length  and  no  fringes  appear  at  Lz .  We  were  able  to  accomplish  this 
with  arc  lamp  sources  (after  considerable  struggle).  The  diffuser  was  a  device 
to  preclude  lnadvertant  "cheating"  by  passing  spatial  Information  from  the  mask 
plane  (at  Mx)  to  the  receiver.  The  receiver  Is  a  simple  Michelson  Interferometer. 
Fringes  obtained  at  the  output  (e.g.  Figure  3.4)  were  exactly  those  expected  for 
the  source  and  geometry  as  calculated  from  the  formula  for  Sm(x)  In  Sec.  3.3. 

Our  first  attempt  to  modulate  Sm(x)  was  a  total  failure.  We  then  saw  that  our 
formulas  and  (in  retrospect)  common  sense  predicted  that  failure.  Because  these 
results  are  obvious  to  us  only  in  retrospect .  we  believe  It  to  be  Important  to 
review  them  here.  We  attempted  to  modulate  the  envelope  (visibility  curve  at  the 
output  by  blocking  certain  parts  of  the  visibility  curve  with  a  mask  In  the  plane  of 
M4.  The  result  of  this  action  decreased  the  total  transmission  but  had  no  effect 
on  the  output  visibility  curve.  It  Is  now  clear  why  this  Is  true.  Our  masks  blocked 
many  fringes,  so  It  obviously  blocked  just  about  as  much  of  one  wavelength  as  of 
any  other  .  Only  masks  containing  detail  of  about  the  size  of  a  fringe  avoid  this. 

In  order  to  have  a  significant  effect  on  a  system  using  light  of  spectral  width  Act  , 
the  length  of  the  mask  must  encompass  a  path  mismatch  well  In  excess  of  the  value 
l/Acr  needed  to  barely  resolve  that  width.  Thus  we  need  a  mask  filling  M4  with  a 
lot  of  fringe-size  detail  In  It  to  give  good  modulation.  We  would  need  to  use  stan¬ 
dard  Integrated  circuit  photolithography  methods  to  make  the  masks  described 
analytically  in  Sec.  3.3.  Therefore,  we  omitted  the  actual  demonstration  as  well 
beyond  our  limited  funding  level. 

3 . 5  CONCLUSION  ON  THE  NEW  METHOD 

The  new  spectral  domain  spread  spectrum  method  has  been  shown  to  work,  although 
In  a  slightly  more  complicated  way  than  we  had  anticipated. 

The  receiver  compi  lses 

•  a  double  Interferometer, 

•  N  predesigned  spatial  masks  l  aligned  normal  to  the  x  direction  In 
the  spatial  Fourier  transform  plane, 

«  N  acousto-optic  modulators  (one  for  each  of  the  masks) . 

Each  modulator  transmits  its  own  Independent  message.  The  outgoing  beam  is 
uniphase  and  spatially  uniform. 
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Figure  3.4  The  "White  light"  Spectral  Fourier  Transform  As 
Observed  In  The  Spectral  Filter  Plane 
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The  receiver  comprises 


.  N  detectors  (one  at  each  N  preselected  location  In  the  Fourier  transform 
of  the  received  spectrum. 

Each  detector  reads  the  message  sent  by  the  corresponding  modulator. 

Many  problems  remain,  but  the  basic  Idea  Is  now  understood  and  known  to  be  both 

feasible  and  promising  (for  jam  resistance,  security,  etc.). 

4  CONCLUSIONS 

This  study  has  led  to  certain  definite  conclusions.  In  particular 

(1)  Optics  can.  In  principle,  be  used  for  synchronization  of  very  long  signals  even 
though  using  such  signals  may  not  be  very  meaningful. 

(2)  For  signals  longer  than  about  1000  bits  but  not  so  long  as  to  Invoke  the  limita¬ 
tion  just  mentioned,  optical  synchronization  Is  very  appropriate. 

(3)  It  Is  feasible  to  do  all-optical  decoding  of  PRC  SSC  signals  while  tracking 
Doppler,  fading,  and  synchronization,  thus  achieving  significant  Improvements 
In  signal-to-noise  ratio. 

(4)  Acousto  optic  tunable  filters  can  operate  at  resolving  powers  orders  of 
magnitude  higher  than  previously  supposed.  This  has  important  applications 
to  Doppler  tracking  and  to  the  decoding  of  frequency  hopping  systems . 

(5)  Matched  filtering  Is  a  special  case  of  a  far  more  powerful  technique  which  we 
call  "generalized  matched  filtering". 

(6)  Matched  filtering  is  a  useful  tool  for  burst-error-correcting  decoding. 

(7)  Fourier  transform  spectral  multichannel  transceivers  have  been  designed  to 
°;ive  extremely  great  security  of  comml cation. 

The  net  effect  of  these  conclusions  Is  to  emphasize  the  truth  of  our  Initial  assump¬ 
tion  that  optics  can  be  an  Important  part  of  spread  spectrum  signal  decoding. 
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